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INTRODUCTION 


u 

This is the September 1995 Semi-Annual report for -NA45bW^3023-(4>w.R4*Rroj 
Studies of Extm-Solar Oort Clouds and the Kuiper Disk, (S.A. Stern, PI). 

We are conducting research designed to enhance our understanding of the evolution and 
detectability of comet clouds and disks. This area holds promise for also improving our under- 
standing of outer solar system formation, the bombardment history of the planets, the transport 
of volatiles and organics from the outer solar system to the inner planets, and to the ultimate fate 
of comet clouds around the Sun and other stars. According to “standard” theory, both the Kuiper 
Disk and the Oort Cloud are (at least in part) natural products of the planetary accumulation stage 
of solar system formation. One expects such assemblages to be a common attribute of other solar 
systems. Therefore, searches for comet disks and clouds orbiting other stars offer a new method 
for inferring the presence of planetary systems. 

This project consists of two major efforts: (1) observational work to predict and search for 
the signatures of Oort Clouds and comet disks around other stars; and (2) modelling studies of ' 
the formation and evolution of the Kuiper Disk (KD) and similar assemblages that may reside 
around other stars, including (3 Pic. These efforts are referred to as Task 1 and 2. The main 
collaborators with PI Stern in Task 1 are Cols Drs. David Weintraub (Vanderbilt U.) and Mike :j 
Shull (U. Colorado). The main collaborator in Task 2 is Col Dr. Glen Stewart (U. Colorado). j 

RECENT PROGRESS 

Task 1: Observational Studies of Comet Disks and Clouds 

Last year we completed 2 new observing runs at the IRAM and SEST submm telescopes to 
study one of the best IRAS IR-excess comet cloud candidates, a PsA (Fomalhaut), r Ceti, (3 Pic, 
and HD98800. Our HD98800 measurements were judged of immediate interest and published in 
an IAU Circular (IAUC 6003; June 1994). 

Our observing work will continue in early 1996 with multi-channel bolometer observations 
^ of IR excess stars at IRAM and (if approved) a CSO run to survey ~ 15 pulsars for evidence of 

cold dust emission indicative of planetary formation. 

In parallel with our observing program, two papers on theoretical subjects related to comet 
cloud and planetary detection have been completed, submitted, and published. In the first, PI 
Stern developed a new method for directly detecting planets. This method is based on the fact 
that giant impacts, like the one forming the Earth-Moon system, cause the target planet to become 
IR-luminous for several hundred to several thousand years. The paper demonstrates (i) that such 
events will be detectable from the twin-Keck interferometer, and (ii) that in a young star cluster 
like Orion, between one and a few such “hot planets” should be extant at any given time. In the 
second paper, PI Stern and collaborator Shull have shown that the lack of a galactic population of 
Gamma Ray Bursters (GRBs) constrains the fraction of stars in the galaxy that have Oort Clouds 
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to < 20%. These papers have now appeared in print. 

Task 2: Collisional Modelling of the Kuiper Disk 

We have now completed the first model of collision rates of the Kuiper Disk. With this 
model we explored the rate of collisions among bodies in the present-day Kuiper Disk as a function 
of the total mass and population size structure of the Disk. We find that collisional evolution is an 
important evolutionary process in the Disk as a whole, and indeed, that it is likely the dominant 
evolutionary process beyond « 42 AU, where dynamical instability timescales exceed the age of 
the solar system. Two key findings we report from this modelling work are: (i) That unless the 
Disk’s population structure is sharply truncated for radii smaller than ~l-2 km, collisions between 
comets and smaller debris are occurring so frequently in the Disk, and with high enough velocities, 
that the small body (i.e., km-class object) population in the disk has probably developed into a 
collisional cascade, thereby implying that the Kuiper Disk comets may not all be primordial, and 
(ii) that the rate of collisions of smaller bodies with larger 100 < R < 400 km objects (like 1992QBi 
and its cohorts) is so low that there appears to be a dilemma in explaining how QBjs could have 
grown by binary accretion in the disk as we know it. Given these findings, it appears that either 
the present-day paradigm for the formation of Kuiper Disk is failed in some fundamental respect, 
or that the present-day disk is no longer representative of the ancient structure from which it 
evolved. In particular, it appears that the 30-50 AU region of the Kuiper Disk has very likely 
experienced a strong decrease in its surface mass density over time. This in turn suggests the 
intriguing possibility that the present-day Kuiper Disk evolved through a more erosional stage 
reminiscent of the disks around the A-stars /? Pictorus, a PsA, and a Lyr. These results were 
published in The Astronomical Journal this year. 

We have also used this model and a second code to estimate the detectability of IR emission 
from debris created by collisions. We found that eccentricities in the Kuiper Disk are high enough 
to promote erosion on virtually all objects up to ~ 30 km, independent of their impact strength. 
Larger objects, such as the 50-170 km radius “QBi” population, will suffer net erosion if their 
orbital eccentricity is greater than « 0.05 (« 0.1) if they are structurally weak (strong). The 
model predicts a net collisional erosion rate from all objects out to 50 AU ranging from 3 x 10 16 to 
10 19 gyr, -1 depending on the mass, population structure, and mechanical properties of the objects 
in the Disk. We find two kinds of collisional signatures that this debris should generate. First, there 
should be a relatively smooth, quasi-steady-state, longitudinally isotropic, far IR (i.e, ~ 60 fi m 
peak) emission near the ecliptic in the solar system’s invariable plane ecliptic, caused by debris 
created by the ensemble of ancient collisions. The predicted optical depth of this emission could be 
as low as 7 X 10 -8 , but is most likely between 3 X 10~ 7 and 5 X 10 -6 . We find that this signature 
was most likely below IRAS detection limits, but that it should be detectable by both ISO and 
SIRTF. Second, very recent impacts in the disk should produce short-lived, discrete clouds with 
significantly enhanced, localized IR emission signatures superimposed on the smooth, invariable 
plane emission. These discrete clouds should have angular diameters up to 0.2 deg, and annual 
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parallaxes up to 2.6 deg. Individual expanding clouds (or trails) should show significant temporal 
evolution over timescales of a few years. As few as zero or as many as several 10 2 such clouds 
may be detectable in a complete ecliptic survey at ISO’s sensitivity, depending on the population 
structure of the Kuiper Disk. This work was recently accepted for publication in Astronomy & 
Astrophysics. This paper is attached in preprint form. 

These two papers were also accompanied by an invited review, submitted to the Planetary 
Ices book, summarizing the present state of knowledge about the Kuiper Disk and Pluto. Addi- 
tionally, PI Stern gave three invited talks summarizing the collisional modelling results obtained 
under the Origins program. A list of these invited talks is attached. 

Finally, in July 1995, we organized and sponsored a 2-day workshop on collisions in the 
Kuiper Disk. This workshop was attended by D. Davis (PSI), P. Farinella (Italy), R. Canup 
(U. Colorado), M. Festou (France), J. Colwell (U. Colorado), H. Levison (SwRI), and PI Stern 
(SwRI). The proceedings of this workshop were informally published and the distributed among 
the participants. A copy was also sent to Origins program scientist Trish Rogers. 
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Abstract 


We explore the rate of collisionally-produced mass injection in the present-day Kuiper Disk 
as a function of the total mass and population size structure of the disk. Our objective is 
to estimate the detectability of IR emission from debris created by collisions. We find that 
eccentricities in the Kuiper Disk are high enough to promote erosion on virtually all objects 
up to ~ 30 km, independent of their impact strength. Larger objects, such as the 50-170 km 
radius “QBi” population, will suffer net erosion if their orbital eccentricity is greater than 
» 0.05 (« 0.1) if they are structurally weak (strong). Our model predicts a net collisional 
erosion rate from all objects out to 50 AU ranging from 3 x 10 16 to 10 19 g yr -1 , depending 
on the mass, population structure, and mechanical properties of the objects in the Disk. We 
find two kinds of collisional signatures that this debris should generate. First, there should 
be a relatively smooth, quasi-steady-state, longitudinally isotropic, far IR (i.e, ~ 60fJ.m 
peak) emission near the ecliptic in the solar system’s invariable plane ecliptic, caused by 
debris created by the ensemble of ancient collisions. The predicted optical depth of this 
emission could be as low as 7 x 10~ 8 , but is most likely between.3 x 10 -7 and 5 x 10~ 6 . We 
find that this signature was most likely below IRAS detection limits, but that it should be 
detectable by both ISO and SIRTF. Second, very recent impacts in the disk should produce 
short-lived, discrete clouds with significantly enhanced, localized IR emission signatures 
superimposed on the smooth, invariable plane emission. These discrete clouds should have 
angular diameters up to 0.2 deg, and annual parallaxes up to 2.6 degrees. Individual 
expanding clouds (or trails) should show significant temporal evolution over timescales of 
a few years. As few as zero or as many as several 10 2 such clouds may be detectable in a 
complete ecliptic survey at ISO’s sensitivity, depending on the population structure of the 
Kuiper Disk. 
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1. Introduction 


Almost, a half- century ago, Edgeworth (1949) and Kuiper (1951) made prescient pre- 
dictions that the Sun should be surrounded by a. disk-like ensemble of comets and other 
“debris” located beyond the orbit of Neptune. The case for such a primordial reservoir 
was strengthened when it was later pointed out that such a disk could be an efficient 
source region to populate the low-inclination, short-period comets (i.e., the Jupiter Family 
Comets; Fernandez 1980). Convincing dynamical simulations supporting this link between 
the Jupiter Family comets and the Kuiper Disk (KD) region arose when Duncan, Quinn, 
& Tremaine (1988) and later Quinn et al. (1990) showed that a low-inclination source 
region appears to be required for a low-inclination orbit distribution of the Jupiter Family 
Comets. Figure 1 is a schematic diagram revealing the gross architecture of the this disk 
in relationship to the orbits of the five known outer planets. 

The computational capabilities available to Duncan and coworkers in the mid-1980s 
required some important approximations be accepted (for reviews of this work, cf., Weiss- 
man 1993, and Stern 1995a). Although these compromises were criticized by Bailey & 
Stagg (1990), Duncan et al.’s work generated interest in the Kuiper Disk by both mod- 
ellers and observers. Of particular relevance here are the Holman & Wisdom (1993) and 
Levison & Duncan (1993) studies of orbital evolution in the Disk. These workers found 
a time-dependent dynamical erosion of the Disk population inside ~ 42 AU, caused by 
nonlinear perturbations from the giant planets. The dynamical chaos resulting from these 
perturbations is ultimately responsible for the transport of short-period comets from the 
long-lived Kuiper Disk reservoir to planet-crossing orbits where they can be routinely 
detected. Based on the bias-corrected population of Jupiter Family Comets and the dy- 
namical transport efficiency of comets from the Kuiper Disk to the inner planets region, 
Duncan, Levison, & Budd (1995) have since estimated that ~ 10 10 comets orbit, in the 
Disk between 30 and 50 AU from the Sun. 

Observational confirmation of the Kuiper Disk was first achieved with the discovery 
of object 1992QBi by Jewitt & Luu (1993). As of late 1995, no fewer than 30 QBi-like, 
trans-Neptunian objects have been discovered (Jewitt & Luu 1995; Stern 1995a). These icy 
outer solar system bodies are expected to have dark surfaces consisting of an icy matrix 
contaminated by silicates and organics. Assuming a typical (i.e., cometary) geometric 
albedo of 4%, and the absence of coma, the distances and magnitudes of these objects 
indicate they have radii between roughly 50 and 180 km. Based on the detection statistics 
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obtained to date, one ran easily estimate that a. complete ecliptic survey would reveal 
~ 3.5 x in' 1 such bodies orbiting between ~30 and 50 AU. Simple power law extensions 
of this population predict a cometary population (which we define as bodies with radii 
between 1 and 6 km) of ~ 1() 10 , which is similar to the dynamical modelling estimates 
obtained by Duncan et al. (1995) to satisfy the short-period comet flux. Very recently, 
Cochran et al. (1995a., b) have reported Hubble Space Telescope results giving the first 
direct evidence for comets in the Kuipcr Disk. 

The objective of this paper is to examine whether there might be detectable signatures 
of collisions in the Kuiper Disk. An extensive study of collision rate dependencies in 
the Kuiper Disk has been published by this author (Stern 1995). We are aware of an 
unpublished manuscript by Alcock fc Hut (1995) on the detectability of individual flashes 
from comet-comet collisions in the Kuiper Disk, and an abstract by Dasgupta et al. (1994) 
predicting the small grain population of the Kuiper Disk; the latter has recently been 
expanded into a short paper by Backman et al. (1995). The study presented here is more 
extensive than these previous works, and is aimed specifically at predicting signatures the 
ISO and SIRTF spacecraft might detect. 

This paper is organized as follows. In §2, 1 summarize our previously-published model 
for estimating the rate of collisions in the Kuiper Disk, on which many of the results 
presented here are based. In §3, we evaluate whether collisions in the Disk are primarily 
erosional or accretional, showing that present-day collisions are expected to be erosive, 
thereby injecting debris into the Kuiper Disk. In §4, I summarize some key results derived 
from that model. In §5, I estimate the quasi-steady-state debris mass, debris opacity, and 
IR brightness expected from the disk. In §6, I discuss the characteristics and number of 
discrete, expanding debris “clouds” produced by individual collisions during the past few 
years. In both §5 and §6, comparisons to IRAS, ISO, and SIRTF detection capabilities 
are made. To close, §7 summarizes the various results and describes how the Kuiper Disk 
debris detection can make important contributions to our understanding of the mass, size 
structure, and mean random eccentricity of objects in the Disk. 


2. The Collision Rate Model 

As described more fully in Stern (19951); hereafter, Paper I), our model for estimating 
collision rates in the Kuiper Disk begins by defining the Disk in terms of a power law 


3 



exponent., o, on the size distribution of objects ill the Disk. Defined in this way, the 
number of objects dN(r) between radius r and r + dr is given by 

dN(r) = Nor" dr, (1) 

where No is a normalization constant set by the estimated number of QBi objects. We 
treat this size distribution as a monotonically decreasing series of radius bins. In what 
follows, we choose bodies in each successive radius bin r to be a factor of 1.6 times larger 
in size (and equivalently, 4 times higher in mass), ranging from r = 0.1 km to r = 162 km. 

We also define a. power law exponent (3 on the radial distribution of heliocentric surface 
mass density £(/i) in the disk, so that 


T,{R) = T,oR p , 


( 2 ) 


where So is the normalization constant. Because there is presently no information on 
the way in which ensemble-averaged inclinations ((z)) and eccentricities ((e)) vary in the 
Kuiper Disk, we adopt a disk-wide (i) and (e) for each run, and vary these quantities 
from run to run as free parameters to explore how sensitive the model results are to these 
variables. 

Once the disk is defined as described above, we bin the disk into a series of concentric 
tori (at successive semi-major axes a) that are 1 AU in width. For each size bin/heliocentric 
bin pair, the model computes the collision rate that a target of given size experiences from 
potential impactors in bins of equal or smaller size. This is of course a function of (ae), 
since (e) controls both the internal velocity dispersion in the disk, and (ae) controls the 
degree of heliocentric bin crossing. 

To compute collision rates we adopt a particle in a box formalism. This approach 
states that the instantaneous collision rate c of target bodies with semi-major axis a, 
eccentricity e, and radius r x being struck by impactors of radius r y , is just 

c(r x ,r y ,a,e) = na g v, (3) 

where n is the local space density of impactor bodies, v is the local average crossing 
velocity of the target body against the impactor population at distance R, and <r g is 
the gravitational-focusing corrected collision cross section of the impactor+target pair. 
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Gravitational focusing is an important correction for targets in the QB| size range and 
larger, particularly in the case of very low eccentricities (e.g., (e) < 1() -2 ) . The orbit- 
averaged collision rate c{r x ,r y ,u,e) can he written to show all its explicit dependencies in 
the model as 

<i(l+(e)) 

^ ^ /(a, (e), /?.) n.( ? y, /?.) UjpyCa, ^e) , (f) ,/?) cTylr'j; , J’y, , v e >c(r+j/) )i 

R = a( l-(e)) 

(4) 

where the term / represents the ratio of the time the target body in an orbit defined 
by (a, (e), (i)) spends in each torus it crosses during its orbit ( T(a , (e),7?)), to the target 
body’s orbital period, « / • 

To compute T(a, (e), 71) we solve the central-field, Kepler time of flight equation explic- 
itly for every (a, (e)) pair in the run’s parameter space. The number density of impactors 
n(r y ,/2) in the torus centered at distance R is computed from the mass of the Disk, the 
Disk’s wedge angle (z), its heliocentric surface density power-law, and the Disk’s population 
size structure power law. The collision timescale r(r x , r y , a, e) is just c" 1 (r x , r y , a, e). 

The code which implements this model produces an array containing collision rates 
c(r x , r y , a, e ) throughout the specified disk, where the free parameters defining the disk are 
the total number of QBi’s, o, /?, and (e). We assume { i ) = |(e). 


3. Are Kuiper Disk Collisions Predominately Accretional or Erosional Today? 

Before continuing, it is relevant to ask whether collisions in the Kuiper Disk are ener- 
getic enough to primarily promote erosion, or alternatively, sufficiently gentle to promote 
growth. 

Whether a collision between an impactor and a target results in growth or'erosion 
depends primarily on the energy of the impact and the mass and strength of the target. In 
the Kuiper Disk, the typical approach velocity of two objects at a distance large compared 
to the Hill sphere of the target can be reasonably-well approximated by 

Voo = \[2v k ((e 2 ) + (i 2 )) 5 , (5) 

where vk is the local Keplerian velocity. For (z) = |(e) we have 

Voo = >/3 (e 2 )5vK. (6) 
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The energy at. impact is therefore given by 


> A tw imp 




(7) 


where // is the reduced mass and v esc is the escape velocity of the two colliding bodies 
measured at the radius of impact. The critical velocity for erosion to occur is given 
approximately by the requirement that the specific impact energy must exceed the total 
energy required to break up the surface and then disperse the ejecta out of the gravitational 
well of the impactor/ target collision pair. The impact energy, Ej mp , as given by Eqn (7), 
must equal or exceed these energy sinks. The critical condition for the target to lose mass 
is that the mass of the ejecta exceeds the mass of the impactor. Therefore, if the impactor 
mass is small compared to the target, then one requires 



> 


/v 



1 

+ 2 U 


2 

esc 



( 8 ) 


where v s represents the velocity required to mechanically shatter the target surface, v esc 
represents the velocity required to disperse the debris to infinity, and k is a factor that 
takes into account energy partitioned into heat, sublimation, hydrodynamic effects, and 
other factors. In what follows we assume K = 8 or 10 (cf., Davis et al. 1989). Now we take 
the specific energy for mechanical breakup of the target to be 


v c 


= *, 


(9) 


where s is the specific strength of the target material at zero compression. Of course, 

-/ 


I2GM, 


Vesc 


( 10 ) 


where G is the universal gravitational constant, Mt is the combined mass of the target and 
impactor, and r t is the combined radii of these two objects. Therefore, from Eqns (7)-(10) 
one can derive the condition which we must solve for e*, the critical eccentricity at which 
impact energies are high enough to promote net erosion: 


Vesc + 2\/(3)e*Vtf Vesc + 3(e*) 2 v^ - KV? - kv\ \ = 0 , 


( 11 ) 
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Notice e* is a function of several target parameters, including the strength, size, and mass. 

Table 1 solves Equation (11) solutions for the critical erosion eccentricity e*, both for 
impacts onto strong (e.g., rock/ice) targets (p = 2 g cm -3 and s = 3 x 1() G erg g -1 ), and 
relatively weak (e.g., snow-like) targets (p — 0.5 g cm -3 and s = 3 x 10' 1 erg g _1 ; these may 
be more applicable to comets in the Disk, cf., Asphaug & Benz 1994); results are shown 
at heliocentric distances of 35 and 60 AU. Following the results discussed in Fujiwara et 
al. (1989), we assume v e j = 0.2i >, mp . 

The results presented in Table 1 can be summarized as follows. For targets with radii 
smaller than ^30 km, e* < 0.01, independent of heliocentric distance and the strength of 
the target surface (within our parameter space). For reference, at 40 AU an (e) = 0.01 
corresponds to a typical encounter velocity at infinity of 82 m/s. Since all objects with 
measured eccentricities in the Disk have eccentricities in excess of 0.01 ( Jewitt & Luu 1995; 
H.F. Levison, pers. comm. 1995), we conclude that most objects of radius a few tens of 
kilometers or less in the Disk are likely to be be eroded when they suffer a collision. These 
results remain valid even if the characteristic ejecta velocity is as low as 5% of the impact 
speeds instead of the 20% used in Table 1. 

For QBi-sized targets with radii near 100 km, we find e* > 0.02 0.03 are required for 

net erosion if they are weak, and e* > 0.05 0.06 is required if they axe strong. Similarly, 

for QBi-like objects with R=170 km, which is comparable to the largest-discovered objects 
in the Disk to date, e* > 0.04 — 0.05 is required if the objects are weak and e* > 0.09 — 
0.10 is required to result in net erosion if the objects are strong. Interestingly, it is also 
worthwhile to note that if the characteristic ejecta velocity is as low as a few percent of 
the impact speeds, then e* will rise dramatically and the QBi population will be in an 
accretional regime, even for eccentricities as high as 0. 2-0.5. We thus find that some QBi’s 
should be undergoing erosion, while others may be in an accretional regime, depending on 
their eccentricity and strength. However, until much better eccentricity statistics become 
available, it is not possible to determine if the QBi population as a whole is gaining or 
losing mass. All we can say is that the range of detected eccentricities span the range of 
e* ’s, creating a complex situation. 

In summary, small objects like comets are being eroded when collisions occur on them 
in the Kuiper Disk, but whether larger objects like the QBi’s suffer erosion or accretion 
depends on their surface strength and the relative eccentricity of the impactor/target pair. 
Regardless of the net fate of the QBi population, collisions inject collisions! debris into 
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the present, -clay Disk. 


4. Collision Rates in the Present-Day Kuiper Disk 

In Paper I, an extensive set of collision rate calculations were presented. Here, we 
briefly summarize some results obtained from those calculations as a prelude to the calcu- 
lation of observable parameters arising from those collisions. 

In the work presented in Paper I and what follows below, it is assumed that the Kuiper 
Disk has an inner radius of 35 AU and an outer radius of 70 AU. We let eccentricity range 
as a free parameter from 2 x 10“ 4 to 2 x 10“ 1 , which extends over a range of eccentricities 
considerably greater than actually seen in the Kuiper Disk objects detected to date. As 
noted before, we assume ( i ) = |{e). 

Four cases defining the heliocentric mass dependence and size distribution of objects 
in the Disk will be presented here. These four cases represent the combinations of two size 
distributions of objects (cf., a in Eqn. 1) and two heliocentric mass distributions (cf., /? in 
Eqn. 2). 

Our favored size distribution, which we call the nominal (NOM) case, connects the 
observationally-estimated ~ 3.5 x 10 4 QBi-sized objects (Jewitt & Luu 1995) inside 50 AU 
with the modelling-derived estimates of ~ 10 10 comets (Duncan et al. 1995) in a single 
power law with a = —11/3. Our second case assumes a = —4, which gives a Constant 
Mass (CM) in every logarithmic size bin. Relative to the NOM case which produces ~ 10 10 
for 35,000 QBi’s (100 km in radius or larger), the CM case produces «5x 10 10 comets. 

For the heliocentric distribution of mass in the Disk, the two cases we ran are defined 
as follows. One case assumes a Constant Mass per Heliocentric Bin (CMHB; /? = —1), 
which corresponds to a surface mass density that declines with heliocentric distance as 
R _1 . The second, and more realistic case, assumes a Declining Mass per Heliocentric Bin 
(DMHB; /? = —2), corresponding to a surface mass density falling like R~ 2 . These two 
cases bracket a realistic range of parameter space (Lissauer 1987). 

Table 2 summarizes these four run cases (and a fifth case we shall discuss later). We 
now present some key results. Figure 2 depicts collision timescale results obtained at 40 
and 60 AU for the NOM/DMHB case. The upper panels show the collision timescale for 
(e) ~ 10 -2 . The lower panels show the collision timescale for (e) ~ 10 -1 . These values 
of { e ) are representative of measured eccentricities in the Kuiper Disk (cf., Stern 1995a; 


8 



.1 cwit.l. & Luu 1995). Figure 3 depicts the same kind of data for the CM/DMHB ease. The 
reader is referred to Paper I to see the collision calculations for the CMHB cases. 

As discussed in detail in Paper I, the key results obtained from the. collision model 
are: (i) That unless the Disk’s population structure is sharply truncated below ~1 km, 
comet-comet collisions are occuring so frequently, and with high enough velocities, that the 
primordial cometary population would be significantly depicted over the age of the solar 
system, (ii) That the rate of collisions of smaller bodies with QBj-size (i.e., 100 < R < 600 
km) objects is so low today that there appears to be a dilemma in explaining how QBi’s 
could have grown by binary accretion in the Disk as we know it. And (iii) that collisional 
evolution is important not only at the inner edge of the Disk near 40 AU, but also much 
farther out, at 60 AU, if the Disk extends that far. These findings suggest that collisional 
evolution plays an important role throughout the Ivuiper Disk, and that collisions play a 
dominant role beginning near- 42 AU, where dynamical timescales exceed the age of the 
solar system. Because collisions in the Disk are modifying the size distribution, these 
findings also argue that the population structure of the present-day Disk appears to be 
in disequilibrium. We refer the reader to Paper I for a more complete discussion of these 
aspects of collisional evolution in the Kuiper Disk, and now turn to the primary subject 
of this report — the observable signatures created by collisions in the Disk. 


5. The Smooth Ecliptic Signature 

As described above, we have found that typical collisions taking place today among 
small bodies in the Kuiper Disk axe erosive. It is well known from studies of the asteroid 
belt that such collisions produce a power-law distribution of debris ranging from a largest 
fragment of radius scale a significant fraction of the target size, down to fine dust (e.g., 
Holesapple 1993). Our interest is in computing the optical depth and thermal emission 
brightness of the debris. In what follows we study the collisionally-generated debris tail of 
objects down to sizes as small as 3 fim. 

We begin by first estimating the steady-state amount of debris in the Kuiper Disk, 
which we compute according to 

M d (a) k M(a) (12) 

where M is the time-averaged rate of debris production in the disk, and T[ oss is the loss 
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timescale for such debris. The relevant, loss processes are dynamical transport owing to 
perturbations by the giant, planets, dynamical loss to Poynting- Roberson (P-R.) drag, and 
physical loss to collisions which destroy grains. Radiation pressure is not as effective as P-R 
drag for particles several microns in size and larger (e.g., Burns et al. 1979), so we neglect 
its effects in the scaling calculations that follow. We also neglect particle sublimation, 
because H^O/silicate grains and larger particles in our size range beyond ~15 AU have 
sublimation lifetimes longer than the age of the solar system. 

In our model we take rrJ, = + T- 1 + T d -\ where T<f yn is the dynamical loss 

time for objects suffering perturbations by the giant planets, T pr is the Poynting- Roberson 
timescale for dust and small particles, and T^c is the collisional fragmentation timescale 
for self-collisions among debris and debris collisions with larger objects, where “debris” 
particles are all objects smaller than 0.1 km, down to a 3 gm lower size limit. The Td yn 
data were kindly supplied by H. Levison. 

We compute M(a) by summing the ejected mass over all target-impactor collisions. 
We determine the mass ejected by a given impactor striking a target at velocity gvnp by 
balancing the energy budget of the impact against energy losses to mechanical shattering, 
and escape from the collision-pair’s potential well. We determine the excess velocity of 
the escaping mass following previously published results (cf., Fujiwara et al. 1989); the 
characteristic ejecta velocities we employ are in the range 0.05v tmp to 0.20v, mp . We assume 
that only 10% of the energy of the impact is partitioned into the kinetic energy of the debris. 
This approach is similar to that described in Eqns (7)-(10), but instead of solving for e*, 
we solve for the ejected mass m e j at a specified impact velocity defined by the assumed 
eccentricity of the population and the mutual escape speed of the target+impactor pair. 
The code caps the maximum ejected mass to be the combined mass of the impactor -ftaxget 
pair. The mass ejected in this way is assumed to follow a standard debris fragmentation 
power law from a largest fragment down to fine dust a few microns in size (e.g., Holsapple 
1993). 

We compute the Poynting-Roberson (P-R) drag timescale for a dust particle according 
to 

T pr = 7 x 10 6 r p p p («2 - al)( 1 - e) years, (13) 

(e.g., Gustafson 1994), where r p is the particle radius in cm, p p is the particle density in 
g cm -3 , a 1 is the particle’s initial semi-major axis in AU, and a 2 is the particle’s final 
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radius in AU; in wlia.t. follows we take a-> = 30 AU, where wo assume Neptune’s dynamical 
influence spatially truncates the Kuiper Disk debris distribution. Equation (13) implies 
a. characteristic P-R loss t imescale of 5 x 10° years for a 10-micron ice particle with unit 
density at 40 AU; at GO AU the timescale for such a particle is near 2 x 10 7 years. 

Now consider the debris collision timescale, Td c ■ This is the timescale for debris 
particles to be lost to collisions with one another or by striking a large Kuiper Disk object. 
We compute 2^ c using the same particle in a box formalism as described in §3, with the 
number density of debris particles calculated from the M(a) and the assumption that the 
population N(a) /v d 

During the time interval between their creation and their loss, particles created by 
collisions will experience azimuthal smoothing due to regression of their arguments of 
perihelion, regression of their orbit nodes, as well as Keplerian shear resulting from post- 
collision velocity dispersion, radiation pressure, and differential P-R drag. As it turns 
out, the Kepler shear time for orbits 1 AU apart to separate by 180 deg is very short 
(~ 10 4 years) compared to both the nodal regression timescales (each ~ 10 7 years) and 
the PR drag timescale (~ 6 x 10 5 years). As a result, we predict the development of a 
longitudinally-smoothed, quasi-steady-state distribution of collisional debris to develop. 

M calculations performed using the model run results shown in Figures 2 and 3 predict 
time-averaged debris production rates in the range w3x 10 16 to 10 19 g yr' 1 in the Kuiper 
Disk, depending on the characteristic ejecta speed relative to the impact velocity. Given 
this rate of production and the T[ oss timescales discussed above, one sees that the Kuiper 
Disk should contain many comet masses of fine debris at any time. 

Once M(a) and Ti oss {a) are calculated for a given r un case as a function of debris size, 
we can estimate the disk’s quasi-steady state geometric radial optical depth at a specified 
wavelength A in debris of radius r p at heliocentric distance a as 

(14) 

where A d is the integrated area of all the collision products in the heliocentric bin at 
distance a, and A (0( is the total projected area of the Kuiper Disk on the sky, which we 
take as 47ra 2 sin i. A<* is computed according to 

Ad{a,\,r p ) = Ad(a ) r p )7rr 2 T7(A,r p ), (15) 


r KVp) = A igAlp) , 

•Atot(Q) 
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where N^{a,r p ) is the number of debris particles of size r } , at distance a and i/(A,r r ) is the 
optical correction factor for particle cross sections used when r p < A. We adopt i]{ A ) = 1 
unless r p < A, in which case we take 

>/(A)=(^) i0 . (16) 


To obtain a total optical depth r( A), we numerically integrate r(a, A) over r p and a, with 
Nd(a, r p ) provided from the model described here and in §2 and §3. 

Figures 4a. and 46 present optical depth estimates calculated in this way for the two 
collision runs presented in Figures 2 and 3, with an assumed IR observation wavelength 
of A = 60 /mi. These optical depth estimates are presented as a function of (e), p, and 
s. Panel 4c presents the same calculation for a Kuiper Disk population that is truncated 
below radii of 1 km, as suggested in Paper I (the details of this population axe given as the 
final case described in Table 2). Because the population of colliding bodies is so sharply 
reduced in the panel c case, panel c likely gives a lower limit to the 60 pm optical depth 
of the Kuiper Disk. 

On examining the curves in Figure 4, one sees there axe three regimes in each plot. 
The first (i.e., leftmost) regime, which is described by the straight-line decrease in r(e) for 
(e) < e*, is simply the optical depth of the macroscopic (i.e., 0.1-200 km radius) Kuiper 
Disk population. The optical depth decreases with (e) in this regime, primarily because 
(i) = |(e). As a result, the area of the disk on the sky increases, and the optical depth 
is diluted. The second (i.e., middle) regime in Figure 4 begins when (e) = e*. Once this 
condition obtains, t jumps up, since debris generation introduces a new source of optical 
depth. Notice the location of this jump depends on the strength of the surface of the KD 
objects. After the r jump near e*, a third regime obtains, where r again declines owing 
to the dilution effect of increasing the disk eccentricity, as described above. 

The main conclusion we draw from Figures 4a and 46 is that the expected IR optical 
depth for a population with characteristic 0.01 < (e) <0.1 will likely be in the range 
3 x 10 -7 to 10~ 5 , but could be as low as 7 x 10 -8 . Although we will adopt the range 
3 x 10~ 7 to 3 x 10 -5 as a typical steady-state optical depth in the Kuiper Disk, it is worth 
noting that catastrophic, erosive collisions on large QBj objects may occur every ~ 3 x 10' 
years (Stern 1995b). When such collisions occur, the debris mass and optical depth in the 
Disk can further increase by a factor of 50 or more for a timescale comparable to Ti oss . 
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The predicted range of typical-state Kuiper Disk optical depths are higher than the 
value of 10~ 8 detected in the Solar System's most prominant asteroidal dust bands (Low 
et al. 1984), but are 10-2000 times lower than the IR optical depths of the disk ensembles 
around Vega, and many other well-known lR.-excess A stars (Backman & Paresce 1993). 

Is this debris detectable? Figure 5 presents estimates of the debris brightness for three 
temperatures that span the 40 to GO K temperatures expected of large grains in simple 
radiative equilibrium at distances of 35 to GO AU from the Sun. (We note that small grains 
can become hotter than 60 K if they are inside 40 AU, but their loss timescale to P-R drag 
is so short that they do not contribute much optical depth, and therefore do not much 
affect the assumption of a characteristic emission temperature between 40 and 60 K). Since 
the actual distribution of temperatures depends on several unknowns, including the grain- 
size histogram, grain emissivities and albedo, and the heliocentric distribution of grains, 
we choose to present flux density curves for isothermal ensembles at 40 K, 50 K, and 60 K; 
given the level of uncertainty in grain parameters at this time, we believe this is preferable 
to integrating through some arbitrary heliocentric temperature profile. The calculations 
presented in Figures 5a and 56 assume optical depths of r = 3 x 10~ 7 and r = 3 x 10~ 6 , 
respectively. Superimposed on the debris black-body curves are reference detection limits 
expected obtainable with IRAS, ISO/ISOPHOT, and SIRTF (cf., Cruikshank et al. 1990; 
Werner et al. 1991; Backman et al. 1995 have examined the detectability of Kuiper Disk 
debris by COBE). 

Figure 5 indicates that IRAS could have marginally detected the Kuiper Disk if the 
debris is 50 K or warmer and exhibits r > 10“ 6 near 60 pm. Now notice in Figure 4 that 
if (e) < 10“ 2 , then optical depths as high as 10“ 5 could obtain. If this were the case, then 
IRAS would have easily detected the Disk. Since no clear detection has been reported, 
despite searches (cf., Backman & Paresce 1993), we conclude that IRAS providesjndirect 
evidence that (e) > 10“ 2 . Additionally, however, it is possible that IRAS could have made 
a weak detection of the disk if (e) « 10“ 2 . In this regard, the intriguing report of cool, ~ 50 
K “zodiacal” dust bands detected in IRAS data (Rowan- Robinson et al. 1991) warrants 
further analysis. 

Concerning ISO and SIRFT, calculations, which are in agreement with those subse- 
quently produced by Backman et al. (1995) indicate that both ISO and SIRTF should 
be capable of detecting the Kuiper Disk debris background at any plausible temperature 
(note: if higher temperatures than 60 K obtain for some fraction of the Kuiper Disk debris, 


13 



then the emission brightness will be higher, strengthening this conclusion). Kuiper Disk 
collisional debris should be concentrated in a band perhaps 10-20 degrees in width, cen- 
tered on the ecliptic, (or more accurately, the invariable plane). This emission from smooth 
Kuiper Disk debris should be distinguishable from zodiacal emission through its (i) nar- 
rower latitudinal distribution on the ecliptic, and (ii) ~ 5 times cooler color temperature 
of 40 to 60 K. If this emission signature is detected with ISO or SIRTF, it will provide a 
powerful diagnostic for studies of collision rates, and by inference, both the heliocentric 
mass and population size structure of the Disk. 


6. Discrete Signatures from Recent Kuiper Disk Collisions 

In addition to the smooth distribution of debris and IR emission resulting from col- 
lisional debris in the Kuiper Disk, one also expects there may be small, discrete collision 
“clouds” generated by recent impacts which have not yet faded into the debris background. 
Given the low-inclination nature of most orbits in the Kuiper Disk, such clouds will most 
likely appear as small expanding debris trails, with their major axis preferentially aligned 
along the invariable plane. 

We have constructed a first-order model to evaluate the number of collision clouds 
that should be detectable. To be conservative, we have assumed that for a cloud to 
be detectable, it must persist with an optical depth ten times that of the Kuiper Disk 
background (cf., §5), and with a A = 60 pm size scale larger than an ISO/ISOPHOT 
angular resolution element (0.8 aremin) for at least 1 year. 

When a given target body such as a comet loses mass in a collision, the debris created 
by the impact will locally enhance the optical depth of the Kuiper Disk. The emission 
flux from such an expanding cloud will peak when the cloud’s optical depth, r c i ^reaches 
unity. The optical depth evolution of the debris cloud can be approximated, for heuristic 
purposes, by assuming the cloud is spherically expanding. The time at which such a cloud 
reaches any specified r c / is given by 


nrc) = 



(17) 


where M c t is the debris mass in the cloud, r p and p p are the characteristic debris particle 
size and density, and v e j is the characteristic expansion velocity of the cloud. As a relevant 
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example, the time at which a cloud of mass equal to a. 2.5 Ian radius comet reaches unit 
optical depth is 


'J'(Tcl)^ 81 


Pd 


1 g cm" 


- 1/2 


( ’> V 

MO microns/ 


- 1/2 


ej 


5 x 10 3 cm s 1 


days. (18) 


The radius of the expanding cloud at this time will be 

- 1/2 


Rci = 7 x 10 9 ( !± \ 

g cm J J WO microns/ 


- 1/2 


cm, 


(19) 


which is about the radius of Jupiter. Such a cloud at a heliocentric distance of 45 AU will 
have an apparent diameter of ~4 arcsec on the sky and will exhibit an annual parallax: of 
2.6 degrees. 

Of course, one need not detect Kuiper Disk collision clouds at r c j = 1. In principle, the 
clouds can be detected against the smoothed Kuiper Disk signal whenever their contrast 
relative to the background is sufficient. In what follows I assume a background optical 
depth of 10 -6 and a cloud:background contrast ratio of 10:1, so that we require r c j = 10“ 5 . 
This gives an expansion timescale of: 

r(r c ,)« 12 ( W ., ) 1/2 (— M )~ 1/2 ( n3 — years, (20) 

\1 g cm A ) \10 microns/ \10 3 cm s -1 / 


with a corresponding radius of 


( Pd > 

. -1/2 , 

( r P ) 

\1 g cm“ 3 / 

' V 10 microns/ 


- 1/2 


cm, 


( 21 ) 


or 0.15 AU. At 45 AU, such a cloud would have an apparent diameter of «11 aremin on 
the sky. 

To improve on these crude scaling calculations, a model has been constructed that 
uses the collision rate estimates described in §4 as its input. For every collision, the model 
calculates the amount of mass ejected, after accounting for energy losses to mechanical 
shattering and escape from the gravity well (cf., Eqns (5)— (11)). Using (i) the ejected 
mass, (ii) a debris population ranging upward from 3 gm in size to the largest fragments 
predicted by each collision pair at the specified (e), and (iii) adopting a collision debris 
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power law like N ( t(r v ) ~ r~ 3 5 dr p (Davis e( al. 1989), the model calculates the spreading 
time to reach a critical r c i = £r, where £ is the cloud :sky contrast factor and r is the 
Kuiper Disk background optical depth computed for the same collision rate data. Given 
the collision rate statistics and spreading times to reach r c /, the model then computes the 
number of clouds expected on the sky. The free key parameters in this model are (e), v e j , 
and the target’s s and p. As before, in agreement with past energy partitioning results 
in catastrophic collisions (e.g., Fujiwara et al. 1989), we vary v e j between 5 x 10 -2 and 
2 x 10 _1 of the impact speed. As usual, we calculate results for two kinds of target bodies — 
strong (ice/rock) targets and weak (snow-like) targets. 

Figure 6 shows the results of such a model calculation using the data in Figure 2 as 
an input, and the assumptions of a A = 60pm background optical depth of r = 10~ 6 . 
Panel a shows the number of collision clouds expected from catastrophic collisions on all 
comets between 1 and 6 km in radius assuming v ej - = 0.20v, mp ; panel b shows the same 
for v e j — 0.05v{mp» Notice: owing to the fact that at lower ejecta velocities a given ejecta 
energy creates more ejecta, the number of clouds increases as the ratio of v e j to v, mp 
decreases. Notice also that the number of clouds falls with increasing (e), which is in part 
due to the decreasing cloud lifetime with increasing u, mp (and therefore increasing v e j). 

Recall that all of the QBi’s with known orbits have e > 10 -2 . For such a population 
the data in Figure 6 predict that the Disk population structure which produces the colli- 
sional timescales shown in Figure 2 should exhibit one or less clouds if the characteristic 
v e j ~ 0.2u, mp . However if v e j ~ 0.05, then this Disk population structure should produce 
a few to perhaps 10 to 50 detectable clouds scattered around the ecliptic, depending on 
the strength of the bodies suffering collisions. 

Figure 7 shows the results of the same model run, using the Disk population structure 
used to produce the collision timescales shown in Figure 3 as an input. For this case the 
collision rates are higher, in turn producing an order of magnitude more discrete clouds 
than in the corresponding case shown in Figure 6. As such, this simulation predicts that 
perhaps as many as 5 recently-created, detectable collision clouds may be scattered around 
the ecliptic if v e j ~ 0.20u, mp , and 30 to 300 clouds if v e j ~ 0 . 05 v tmp ! 

In closing it is worth noting that although the discrete cloud population model data 
shown in Figures 6 and 7 assume £ = 10, it may be that clouds can be detected at much 
lower contrasts, and then verified via their annual parallax. Setting the contrast criteria 
to £ = 2 and retaining the conservative assumption that the background optical depth is 
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r = 10 <J , liicrcJiscs I, lie estimate of detectable clouds by a. factor of just over 2. 


7. Summary 

As described in Paper I, collisional evolution is an important process in the Kuiper 
Disk. We have shown that the mean random velocities in the Disk today (as deduced from 
orbital eccentricities of recently discovered QBj objects) strongly suggest that present-day 
collisions in the Disk are erosive for virtually all objects with radii < 100 km, and for those 
100-170 km QBi objects with eccentricities in excess of ~ 0.05 if they are weak or ~ 0.01 
if they are strong. We predict a time-averaged rate of between 3 x 10 16 and 10 19 g yr -1 
of collisional debris, ranging from multi-kilometer blocks to fine dust is being injected into 
the Disk. Further, we have found that the grinding of objects in the Kuiper Disk produces 
a population of small particles that dominates the optical depth and therefore the thermal 
emission from the Kuiper Disk. 

This paper has identified two types of detectable IR signatures resulting from Kuiper 
Disk collisions. First, as described in §5, the Kuiper Disk should produce a longitudinally 
smoothed, quasi-steady-state population of debris in the solar system’s invariable plane. 
We found that the non-detection of a strong signature of this kind by IRAS indicates that 
the optical depth of this dust at 60 //m is probably < 3 x 10~ 6 . This in turn implies a 
collisional environment in which (e) > 0.01 inside at least ~ 50 AU. For the ( e ) range 
of 0.01 to 0.10 indicated by existing orbit statistics, we redict that the Disk out 50 AU 
produces a total optical depth of 3 x 10~ ' — 5 x 10 -6 and an IR brightness of 0.2 to 5 
MJy/sr at wavelengths between 50 and 100 pm. These results are in agreement with the 
findings of less extensive model calculations described by Backman et al. (1995). The 
Kuiper Disk IR signal should be detectable by ISO. However, if the KD population is 
strongly truncated below size scales of ~ 1 km, then collisions will be less frequent and 
we expect optical depths as low as 7 x 10~ 8 may obtain; this value of r corresponds to an 
IR brightness of 0.07 MJy/sr, which is close to, but still probably above ISO’s detection 
limits. In either case the thermal emission from the Kuiper Disk will be dominated by 
color temperatures in the range 40 to 60 K. 

In addition to the smoothed ecliptic signature, recent collisions which significantly 
erode comets and larger bodies should superimpose local optical depth and therefore IR 
brightness fluctuations on the smooth Kuiper Disk IR background described above. These 
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discrete clouds should have brightness temperatures that are much colder than the zodiacal 
dust, and like the I\D ecliptic background. Such clouds should also appear as expanding 
spherical clouds or elliptical trails with their major axis preferrentially oriented along the 
invariable plane. Before they fade below a 10:1 clouckbackground contrast ratio, these 
clouds should grow to angular scales of many arc minutes, and should display an annual 
parallax of 1 to 2 degrees, depending upon their distance. The timescale for a given collision 
cloud to expand and dilute its optical depth below a 10:1 contrast ratio is of order many 
years, so a sample of many clouds would display a range of sizes and local brightnesses. 

The number of such clouds depends strongly on the population of kilometer and sub- 
kilometer sized objects, the mechanical strength of these bodies, and the mean random 
eccentricity of this population. Our estimates of the population of such discrete clouds 
with a brightness 10 times that of the smooth ecliptic signal range from < 1 up to > 10 2 . 
This range of uncertainty is in turn due to uncertainties in the population structure and 
mechanical properties of the Kuiper Disk population. The clear implication of this result 
is that the surface density of detected clouds on the sky can be used as a probe of these 
otherwise difficult to ascertain properties of the Kuiper Disk population. 
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Table 1: Critical Eccentricities (e*) for Erosion 


Target Radius 

35 AU 
(Strong) 

60 AU 
(Strong) 

35 AU 
(Weak) 

60 AU 
(Weak) 

001 km 

7 x 10~ 3 

9 x 10~ 3 

1 x TO" 3 

1 x 10“ 3 

010 km 

6 x 10~ 3 

7 x 10" 3 

2 x 10" 3 

3 x 10“ 3 

100 km 

5 x 10~ 2 

6 x lO" 2 

2 x 10 — 2 

3 x lO" 2 

170 km 

9 x 10~ 2 

1 x 10" 1 

4 x lO" 2 

5 x lO' 2 


Notes: Strong target implies p = 2 g cm 3 and s = 3 x 10 6 erg g -1 ; weak target implies 
p = 0.5 g cm -3 and s = 3 x 10 4 erg g -1 . In both cases we take k = 8 and r; e j = 0.20u eS c; 
see §3 for additional details. 
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Table 2: Collision Run Cases 


Model 

Mdiifc 

30 AU < R <50 AU 

Population 

Type 

Disk 

Type 

Nqb\s 

30 AU < R <50 AU 

N 

1 ' comets 

30 AU < R <50 / 

CSD le4 25 

0.12 M® 

NOM 

DMHB 

27,740 

7 x 10 9 

CSD le4 23 

0.32 M® 

CM 

DMHB 

31,316 

3 x 10 10 

CSD le4 13 

0.42 M® 

CM 

CMHB 

41,162 

5 x 10 10 

CSD le4 15 

0.1G M® 

NOM 

CMHB 

36,461 

9 x 10 9 

CSD 3e4 15 

0.0/' M® 

NOM 

CMHB 

17,950 

4 x 10 9 


Notes: M^k is the integral mass over all size bins. CMHB=Constant Mass per Heliocentric 
Bin (/? = —1); DMHB=Constant Mass per Heliocentric Bin (/? = —2). NOM=Nominal 
mass dependence per size bin (cc = — 11/3); CM=Constant Mass per size bin (or = —4). 
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Figure Captions 


Figure 1: Schematic depiction of the Kuiper Disk and the orbits of the outer planets, 
including Pluto. The clearing between the orbit of Neptune and the inner edge of the 
present-day disk is created by the dynamical perturbations of the giant planets. The 
location of the outer boundary of the K wiper Disk is not well constrained, and may well 
extend much farther than shown. 

Figure 2: Contours of Kuiper Disk collision timescales (in years) as a function of target 
and iinpactor size for the model run with a mass per heliocentric bin that declines like 
R _1 , and a “Nominal” population size structure (cf., Table 2). Notice, for example, that 
large objects are rarely struck by impactors greater than a few kilometers in radius. 

Figure 3: Contours of Kuiper Disk collision timescales (in years) as a function of target 
and impactor size for the model run with a mass per heliocentric bin that declines like R _ 1 , 
and constant mass (CM) per bin population size structure (cf., Table 2). Again notice that 
large objects are rarely struck by impactors greater than a few kilometers in radius. 

Figure 4- Estimated quasi-steady-state radial optical depths of debris in the Kuiper Disk, 
as a function of (e). Panel a corresponds to the collision rate computation shown in Figure 
2; recall that the population size structure used in this case extended down to objects 
with radii of 0.1 km. Panel b corresponds to the similar collision rate case that generated 
Figure 3. In each panel, the bold line corresponds to debris production from strong (i.e., 
ice/rock) targets; the dashed line corresponds to debris production asstuning the targets 
are mechanically weak (i.e., snow-like). Panel c shows a similar calculation, based on 
the final run case described in Table 2, which has the population size structure sharply 
truncated for sizes below 1 km. The calculations in all panels assume v e j = 0.20i>,- mp . 

Figure 5: Estimated thermal emission brightness for a quasi-steady-state debris optical 
depths of t = 3 x 10 -6 (panel a) and r = 3x 10~ 7 (panel b ), assuming 2*6™ =40 K, 
50 K, and 60 K. Superimposed are the IRAS, ISO/ISOPHOT, and SIRTF sky brightness 
sensitivity limits. 

Figure 6: The estimated number of detectable collision clouds for the collision rate data 
presented in Figure 2, assuming a background optical depth of r=10 -6 , and a cloud: sky 
contrast ratio of £ =10:1. The results in panel a assume v e j = 0. 20 Vi m p', the results in 
panel 6 assume v e j = 0.05u, mp . 
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b'igure 7: The estimated number of detectable collision clouds for the collision rate data 
presented in Figure 3, assuming a background optical depth of r=]0~ G , and a cloucksky 
contrast ratio of (, =10:1. In panel a we assume v c j = 0.20 w, in panel b we assume 
v c j =0.05 Vt'f n p. 
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